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We study the combined effect of a conical topological defect and a Coulomb charge impurity on 
the dynamics of Dirac fermions in gapped graphene. Beyond a certain strength of the Coulomb 
charge, quantum instability sets in, which demarcates the boundary between sub and supercritical 
values of the charge. In the subcritical regime, for certain values of the system parameters, the 
allowed boundary conditions in gapped graphene cone can be classified in terms of a single real 
parameter. We show that the observables such as local density of states, scattering phase shifts 
and the bound state spectra are sensitive to the value of this real parameter, which is interesting 
(DJ0| from an empirical point of view. For a supercritical Coulomb charge, we analyze the system with a 

^ . regularized potential as well as with a zigzag boundary condition and find the effect of the sample 

' topology on the observable features of the system. 

I. INTRODUCTION 

D , The dynamics of Dirac fermions in a 2+1 dimensional conical space-time [1[ or in the presence of a cosmic string 
^ P exhibits a variety of rich nonperturbative quantum features. In spite of strong theoretical interest, the quantum 
properties of such 2D fermionic systems are difficult to observe in the laboratory. The experimental fabrication of 
(T*) ■ monolayer graphene in 2004 |3|-|5|[, whose low energy excitations behave like negatively charged fermions satisfying a 
^ \ two dimensional Dirac equation | 6 [ -fl2l , offers new possibilities to study the effect of topological defects in such lower 
. dimensional fermionic systems |13l-l32]. The Dirac type excitations in pristine graphene are gapless. However, various 
' impurities, electron-electron interactions, substrate structures and other short distance effects can violate the sublattice 
symmetry in graphene, leading to a mass gap which has attracted both theoretical [33-43] and experimental [44-[5Q| 
attention. Thus, the gapped graphene system provides a unique template to study the nonperturbative quantum 
\ features of massive Dirac fermions in the presence of a topological defect. 

■ In graphene the Fermi velocity vp ^ lO^m/s, which is approximately 300 times smaller than the velocity of light. 
' Consequently, a relatively small external Coulomb charge impurity Ze ^ 1 leads to strong nonperturbative electric 
field effects in graphene [361, l5ll - [65| . In a gapped graphene system, the external Coulomb charge is said to reach the 
critical value when the system dives into the negative energy continuum [58, 61, 64] and quantum instability sets in. 
Any given external charge in gapped graphene can therefore be classified as either sub or super critical. These two 
^ ■ different regimes are characterized by markedly different behaviour of the observables such as the local density of 
_C3 ; states (LDOS) [58|. 

In this paper we shall study the combined effect of a conical topological defect and an external Coulomb charge 
impurity on the low energy quantum dynamics of quasiparticles in gapped graphene. When a cone is formed from a 
graphene sheet, the topological defect introduced in the system gives rise to some nontrivial holonomies [HI, HI, |22| . 
The boundary conditions associated with the holonomies can be realized by introducing a suitable flux tube, analogous 
to a cosmic string, passing through the origin [1, 2, 66-70]. In our analysis, such a flux tube shall be used to model 
the conical topological defect on the 2D graphene sheet. Let us now consider the effect of an external Coulomb charge 
impurity in such a system, whose strength could be either subcritical or supercritical. For a subcritical Coulomb 
charge impurity in the presence of the flux tube, we shall show that the quantization of the gapped graphene system 
is not unique and an additional parameter is required to fully characterize the boundary conditions at the origin. 
In order to understand the physical meaning of such a boundary condition, recall that the Dirac description in 
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graphene is valid for low energy or long wavelength excitations. On the other hand, the topological defect as well 
as the Coulomb charge can lead to additional short range interactions, which cannot be incorporated as dynamical 
terms in the Dirac equation. The combined effect of the short range interactions due to the topological defect and 
the Coulomb charge impurity can however be encoded in the boundary conditions [ti'-^]. If we further impose the 
natural requirement that the graphene system conserves probability and the time evolution is unitary, then all the 
allowed boundary conditions can be labelled by a single real parameter. This leads to a one parameter quantization 
of the gapped graphene system, analogous to what was obtained for Dirac fermions in 2+1 dimensional gravity with 
a topological defect [1, 2]. For the gapped graphene system, we show that the experimental observables such as the 
LDOS, phase shifts and the bound state energies depend explicitly on the new parameter that labels the allowed 
boundary conditions. 

For a supercritical value of the charge impurity in the presence of the topological defect, we study the system with 
a regularized Coulomb potential and also with a zigzag edge boundary condition. The regularization of the Coulomb 
potential takes care of the finite size of the external charge impurity and allows the bound states of the system to 
dive into the negative energy continuum [HsI. [61^ . The critical charge in gapped graphene cone is renormalized to a 
value higher than that of the gapless case and the value depends on the gap, the cut off parameter, the topology 
of the system and also on the boundary conditions used to obtain the quasibound state spectra in the supercritical 
region. It will be shown that with the increase in gap or cut off parameter the critical charge in presence of zigzag 
edge boundary condition increases more rapidly than in presence of a regularized Coulomb potential. 

This paper is organized as follows. In the next Section we set up the Dirac equation for gapped graphene cone with 
a point charge at the apex. This is followed by the analysis of the spectrum in the subcritical region, where we obtain 
the scattering phase shifts, bound state energies and local density of states (LDOS) and show how these physical 
quantities depend explicitly on the sample topology. Then we discuss the effect of generalized boundary conditions on 
the spectrum. In the next section the analysis of the corresponding spectrum is done in the supercritical region with 
a regularized Coulomb potential and with zigzag edge boundary condition. We end this paper with some discussion 
and outlook. 



II. DIRAC EQUATION FOR A GAPPED GRAPHENE CONE WITH A COULOMB CHARGE 

Graphene has a hexagonal honeycomb lattice structure which is formed by two inter penetrating triangular sub- 
lattices Q-fil, B ^ and B. Assuming only nearest neighbour hopping in graphene and parameterizing the energy 
difference between the sublattices by £ we have the Hamiltonian as [8] 

H = pY. [u\{Ra)Ub{Ra + ^) + uiiRA + ^h)Ua{Ra)] 

+ [U\{RA)UA{RA)-Ui{RA^Ul)UB{RA^U^^^^ (1) 

Here U\ and Ua {U^ and Ub) are the creation and destruction operators for electrons localized on sites A{B) 
respectively. The vectors Ui{i = 1,2,3) connect one A sublattice point to its three neighbouring B sublattice points. 
The hopping parameter /3 is related to the probability amplitude for electron transfer between neighbouring sites [gI- 
[l2[. Though for an ideal single layer graphene e = 0, by breaking the sublattice symmetry a gap can be introduced 
in graphene [s^-IES] and in our following work we shall consider the massive Dirac excitations of a gapped graphene 
cone. 

From the Hamiltonian H we obtain that in gapped graphene the energy eigenvalues are minimum at the six 
vertices of the first Brillouin zone of graphene and they are known as the Dirac points. Among these points, two are 
inequivalent S 0] • We consider them to be situated at the opposite corners of the Brillouin zone and we denote 
their wave vectors by Ki and K2. Thus we can construct the four linearly independent energy eigenstates 0, 0, [TO] 
of the hopping Hamiltonian denoted by \Ki^ A) ,\Ki, B) ^\K2, A) and \K2^B). The pseudospin indices A and B in the 
eigenstates correspond to that sublattice on which the wavefunction has nonzero amplitude and the valley indices Ki 
and K2 in the eigenstates are distinguished by the manner in which the phase of the wavefunction evolves around a 
lattice site having zero amplitude wavefunction [24] (see Fig(T]). It can also be seen from FigUthat the states with 
valley index K2 can be produced by rotating the corresponding states with valley index Ki by 180° [16]. The basis 
is chosen in such a manner that K2 = — Ki. The low energy eigenstates in graphene can be expressed as a linear 
combination of these energy eigenstates multiplied by envelope functions varying slowly on the lattice parameter scale. 

The low-energy properties of the quasiparticle states in graphene near the Dirac point having valley index Ki , can 
be described by the Dirac equation 



i^^ = [-i{aid^ + a2dy) + mas] ^ = 



(2) 
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FIG. 1: The linearly independent energy eigenstates of graphene are shown. Here the solid and empty circles belong to 
sublattice A and B respectively. 1, g and represent the nonzero amplitudes of the wavefunction at the lattice sites where 
they are assigned and the wavefunction has zero amplitude at all the other remaining lattice sites. Here g = exp(i27r/3) and 
g^ — exp(— 'i27r/3). 



where m denotes the Dirac mass generated due to sublattice symmetry breaking, E is the energy eigenvalue and we 
have set = = 1- The Hamiltonian acts on the array of the slowly varying envelope functions ^ = 

The Pauli matrices cri,2,3 act on the pseudospin indices A, B. 

To study the effect of topology on this system, the formation of a graphene cone is considered by introducing local 
defects in the hexagonal lattice structure of graphene [iGj. When a sector is removed from the plane sheet of 
graphene and the two edges of the sector are identified, the frame {cx^Cy} becomes discontinuous across the joining 
line. Therefore we choose a new set of frames {ex'-,ey>} which is rotated with respect to the old frame by an angle 
if = 9 ^ ^ in the counter clockwise direction (see Figl2]). The x' and y' axes are chosen along the eg direction and the 
—Cr direction respectively [15. 16]. 

For this change of reference frame the wave function has to be transformed by exp(i(/:?<j3/2) to keep the form of the 
Hamiltonian the same [Hi [l6|. Thus the conical topology gives rise to nontrivial holonomies for the pseudoparticle 
wavefunctions. When a cone with angle of deficit is formed, where n can take only discrete values 1,2,3, 4, 5, the 
angular boundary condition obeyed by the Dirac spinor as it goes around a closed path is given by 

^{r,e = 2^) = e^2-(i-t)¥^(^^^ = 0). (3) 

Here (r, 0) denotes the polar coordinate of the lattice points. 
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FIG. 2: Formation of cone from plane graphene sheet and rotation of the coordinate system in order to make it continuous. 



When the cone is formed by removing odd number of wedges of angle ^ from the plane graphene sheet and the 
two edges of the removed portion are identified, the adjacent sites on two sides of the identification line belong to the 
same sublattice (see Fig 12]). Thus the bipartite nature of the hexagonal lattice is broken. Also from Fig (T] we can see 
that rotation of a state with valley index Ki by an odd multiple of angle ^ gives the corresponding state with valley 
index K2 with the same sublattice label. Therefore the removal of odd number of wedges of angle ^ gives rise to 
an additional phase shift affecting the valley indices of the wave function in the boundary condition [TII-IitI, [23|. The 
states with valley index K2 will be affected in the same manner as the states with valley index Ki but there will be 
a relative phase difference of 180° between them. Therefore this boundary condition can be described by involving a 
r2 matrix in it where the matrix r2 operate on the valley indices [HI, [H, [76| . When n is even, this off diagonal matrix 
does not play any role and the exponential factor appearing in the boundary condition just gives ±1 depending on 
the value of n. We diagonalize the matrix r2 for all allowed odd values of n. As a result the valley indices of the 
electronic states become mixtures of Ki and K2. Then the angular boundary condition satisfied for all values of n, 
by a branch of electronic states having a fixed Fermi index, is given by [l5|, [13] 

^(r,6> = 27r) = e^27r[±f ao+(i-t)¥]^(^^ 51 = q). (4) 
Here ctq is an identity matrix which acts on the pseudospin indices A. B. = [ ^ ' where is a mixture of 

\ ^B,K' ) 

K\ and i^2- 

The effect of these holonomies can be modelled by introducing a fictitious magnetic flux tube [19] passing through 
the apex of the cone. The magnetic vector potential modifies the boundary condition on a Dirac spinor as 

^(r,6> = 27r) = e^^^ ^'^"^^^(r, 6> = 0). (5) 

Here dl is a line element on the circumference of the cone at a distance r from the apex, i.e. 

dl = ee r(l - -)dO. (6) 

Substituting (|6]) in (|5]) and assuming that the component Ag of the magnetic vector potential is independent of the 
angle 0, we have from Eq.(j4]) 
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Then an external Coulomb charge localized at the apex of the gapped graphene cone can be equivalently described 
by a suitable combination of electric charge and magnetic flux tube ^] . Let us assume that the Coulomb interaction 
strength is a = where Z is the atomic number of the impurity, e is the electronic charge and k is the dielectric 
constant. Replacing the ordinary derivatives in the Hamiltonian by the corresponding covariant derivatives, the Dirac 
equation for the low energy excitations of gapped graphene cone in presence of a Coulomb charge at its apex is given 

by 



Let 



*(-'^)=e( S!'! \ -'''^ (9) 
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where j is half-integer. Substituting (|9]) in (|8]), we obtain that the leading short distance behavior of the wavefunction 
is given by 

^ ^^"^ where 7 = \J - and v = (10) 

V-'- ~ "6^ 

We should note that the angular part of the wave function is different from that of the planer case, due to the choice 
of the reference frame [l6j. From (p!Q|) we can see that when |a| exceeds 7 becomes imaginary. Therefore, the 

eigenstates ^^^(r) and ^^^(r) becomes wildly oscillatory and have no well defined limit as r ^ 0. For massive 
excitations the critical coupling ac corresponds to that value of a for which E = —m. When m = 0, the value of ac 
is equal to the minimum allowed value of u and depending on the magnitude of Dirac mass and boundary conditions 
ac increases gradually from u. It will be shown that the critical coupling for the gapped graphene cone explicitly 
depends on the angle of the cone and also on the product of gap and cutoff parameter. From the expression of u one 
can see that if we consider the expression u = and j = f , then = | for all values of n. Therefore for analyzing 

the effect of topology the angular momentum channel j = f has not been considered. 

Depending on the strength of the external Coulomb charge compared to that of the critical charge of a gapped 
graphene cone with a particular opening angle, the effect of the charge impurity on the cone can be analyzed in two 
separate regions: subcritical and supercritical. 



III. DIRAC EQUATION FOR A GAPPED GRAPHENE CONE WITH A SUBCRITICAL COULOMB 

CHARGE 

In this Section we discuss the bound and scattering state solutions of the Dirac fermions in a gapped graphene cone 
in the presence of an external Coulomb charge impurity. Following [36], consider the ansatz 



^'^{p) = VrnTEe-2p^-2[F{p)^G{p)] (11) 

and 



^'^{p) = V^^^e-^p^-nF{p)-G{p)i (12) 

(i-t) 



where p = 2r]r, r] = ^Jm^ — 7 = \/v^ — ^ v — and total angular momentum j takes all half integer 

values. Using Eqs. (j8j), (|9]), ([11]) and ([Hj) we get 

where Hp denotes the radial Dirac operator. From Eq. (p!3|) we have 



dF ( olE\ ^ f ma\ ^ ^ 
p— + 7 ]^ -[^^ G = 0. (14) 



dp \ T] J \ T] J 
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and 



dG / aE\ ^ 

dp \ T] J 



-u^ 1 F = 0. 

V 



Substituting the expression of G from Eq.([Tl]) in Eq. (p!5]) we have 

pF'' + (1 + 27 - p)F' - (^^-^]F = 0. 



(15) 



(16) 



In order to proceed, for the moment we assume that the wavefunction vanishes at the charge impurity. Solutions 
of Eq. (fT6|) which obey that boundary condition are given by (75| 



Ea 

F(p) = AiM(7---,l + 27,p), 



(17) 



where Ai is a constant. From Eq. (p!4|) we have 



The upper and lower components of the wavefunctions are 



(18) 



and 



^%{p) = AiVm- Ee-^p^-^ 



Ea 

M(7-— -,l + 27,p 



Ea 

M(7-— ,l + 27,p 



fc^Mfl + 7--,l + 27, 



(19) 



7 1 + 7 ,l + 27,p 



Bound states occur when the wavefunctions reduce to polynomials i.e. when 

aE 



7- 



V 



where 



P \ 1,2,3...., 

The corresponding bound state spectra is obtained as 

En 



when > 0, 
when u < 0. 



m sgn(a) 



1 



(20) 



(21) 



(22) 



(23) 



(p+7)2 



Here the energy should be of the same sign (positive or negative) as a because otherwise the value of p will become 
negative and in our range of interest, it is not allowed. 

The solution of Eq. (p!3|) which leads to physical scattering states when |£^| > \m\ is [75| 



F{p) = AiM (^7 - 1 + 27, p^ + A2p-^^M 

From Eq. ([T4|) we have 

(7-f) 



Ea 

-7 ,1 -27,p . 

T] 



(24) 



G{p) 



AiM 1 + 7 



Ea 



A2p-^'*M i 1 



Ea 



,1-27,P . (25) 
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Here the parameter 77 = \Jm^ — is purely imaginary, i.e. r] = — i/c, [36| where k is defined as = \JE^ — vn? . 
Consequently, the variable p also becomes purely imaginary, p = —2ikr. Using the r ^ 00 limit of the scattering 
states the scattering matrix is obtained as 

From Eq. ([26|) it can be seen that the poles of the S matrix determined by (l + 7 — = 1 — p, where p is a nonzero 
positive integer and {l — =0 when > 0, gives back the corresponding bound states as expected. 

A. Generalized boundary conditions 

The Dirac equation discussed in the previous section is valid for low energy or long wavelength excitations. The 
conical defect as well as the Coulomb charge impurity might give rise to short range interactions in the system, which 
cannot be incorporated as dynamical terms in the Dirac equation. However, the combined effect of those short range 
interactions can be taken into account through the choice of suitable boundary conditions. In systems with unitary 
time evolution, there is a well defined prescription due to von Neumann to determine the allowed boundary conditions, 
which is what we shall follow [tiI-It^. 

From Eq.® it can be seen that the Dirac operator H has an angular part and a radial part. The angular part 
operates on a domain Y{0) which is spanned by the antiperiodic functions e*-^^ where j is a half integer and the 
corresponding boundary condition is kept unchanged. The radial Dirac operator Hp, given by Eq. (p!3|) . is symmetric 
in the domain Vq = Cq^{R'^) which consists of infinitely differentiable functions of compact support in the real half 
line and its adjoint operator Hj^ has the same expression as Hp but its domain can be different. Now to determine 
the domain of self-adjointness of the Dirac operator H, consider the equations 

H^^± = ±^M/±, (27) 



where d has the dimension of length. The total number of square integrable, linearly independent solutions of Eq. (j27|) 
gives the deficiency indices for H and they are denoted by n±. For obtaining n±, Eq. (p!3|) is considered with E replaced 
by To understand the significance of these indices we should note that if an operator is self-adjoint, then it is 
expected to have only real eigenvalues. Thus the existence of imaginary eigenvalues in the spectrum is a measure 
of the deviation of an operator from self-adjointness. The non zero deficiency indices serve as the measurement of 
this deviation. Depending on the deficiency indices Hp can be classified in three different ways [ill : (1) When 
n+ = n_ = 0, Hp is essentially self-adjoint in Vo{Hp). (2) When n+ = n_ 7^ 0, Hp is not self-adjoint in Vo{Hp) but 
it can admit self-adjoint extensions. (3) When n+ 7^ n_. Hp cannot have self-adjoint extensions. 
To find the deficiency indices n± let us first consider the following. 



*i±(p) = \h ± y^P''-HF±{p) + G±{p)] (28) 



and 

*i±(/') = ^f^e-ip''-i[F±ip) - G±(p)], (29) 
where p = 2?7ir, rji = 7 = a/^^ — and u = "[j"^- Then we can write 

„ f f±W \ Z "4 + ») -(-+¥) \(FM\_r, ,30, 



From Eq. (|3Q|) we have 

dF^ip) 



P- 



dp 
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and 



Substituting the expression of G±{p) from Eq. (|3T]) in Eq.(j32j) we have 

+ (1 + 27 - p)F!,{p) - (7 T ^) F±{p) = 0. (33) 

We first determine the deficiency subspace characterized by Fj^{p) and given in Eq. (|3Q|) . The required solution 

of Eq.(|33D is 

F+(p) = [/(^7-^,l + 27,p). (34) 
Using the differential recursive relation zU'{a^ 6, z) + all (a, 6, z) — a(l + a — 6)[/(a + 1, 6, z), from Eq. ([3T]) we have 



(iy+2^j V 



In the limit p — !• the functions behave as 



where 



P+ = ^ (38) 

r(-7-^)r(i + 27) ^ ^ 

Q+ (39) 

^ r(7-^)r(i-27) ^ ^ 

^ (^ + ^) r(i-7-^)r(i + 27) ^ ^ 

(^ + ^) r(i + 7-;^)r(i-27) ^ ^ 

are constants depending on the system parameters. From the above relations we find that as p — > 0, 

y \^A^fpdp^ J{aip^^ +a2 + asp-^^)dp, (42) 

J \^B+fpdp^ Jihp^^ +62 + bsp-^^)dp, (43) 

where 0^,6^ (z = 1,2,3) are constants, whose explicit forms are not relevant. As 7 is a real positive quantity in 
the subcritical region, from Eq. (|42]) and (|43l) it can be shown that and are square integrable everywhere 
provided < 7 < ^. Thus n+ = 1 for the parameter range < 7 < ^. 

In a similar way, by analyzing the deficiency subspace characterized by the negative sign, we obtain 

F_=C/(^7+^,l + 27,p), (44) 

G- = -fc^fHz2lMf7fi + ,+ i^,l + 27,pV (45) 



md 
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In addition, in the limit p — > the functions F_ and G_ behave as 

sin7r(l + 27) 

sin7r(l + 27)^ ^ ^' ^ ^ 

where 

P-=P+, Q- = Q+, R-=R+, S-=S+. (48) 

Similar analysis as before shows that n_ = 1 for the parameter range < 7 < | as well. Thus for the gapped 
graphene cone with a charge impurity, n+ = n_ = 1 when < 7 < ^. Therefore, this system admits a one parameter 
family of self-adjoint extensions for < 7 < \. We would now like to find out the spectrum of the system in a range 
of V and the effective subcritical Coulomb potential strength a such that < 7 < |. The deficiency subspaces for the 
radial Dirac operator Hp are spanned by the elements, 

'b± 




*± =( i"" ) = I '1 — - « _ : 1- (49) 



The domain in which the Dirac operator is self-adjoint is given by Vz{Hp) = Vo{Hp) {C(e*t^+ -|- e~*t^_)}, where 
C is an arbitrary complex number and z G R mod 27r. Thus we have a one parameter family of self-adjoint extensions, 
labeled by a real parameter z. For each choice of the parameter we have a domain of self-adjointness of the radial 
Dirac operator defined by Vz{Hp). When p — > an arbitrary element G Vz{Hp) can be written as 




where F- and G- denote the complex conjugates of and respectively. 

Now the spectrum of the system is found out when the boundary conditions are governed by the domain Vz{Hp). 
A solution of the physical eigenvalue problem is written as 



/ VmT^e-ip'^-nF(p) + G(p)) \ , . 

Here F and G satisfy Eqs. ([H]) and ([15]) respectively, and N denotes the normalization. Solutions of Eqs.([TH) and 
([T5]) that are square integrable at infinity are given by 

F = C/(^7-^,l + 27,p), (52) 

G = C/(^l+7-^,l + 27,p). (53) 
Now using Eqs. (|5Q|) and (|5T]) we have in the limit p — > 0, 

sm7r(l + 27) 

G -{R - Sp-^^), (55) 

sin7r(l + 27)^ ^ ^' ^ ^ 



where 



r(-7-^)r(i + 27) r(7-^)r(i-27) 



(56) 



R=—. ^ S = —, . (57) 

r(i-7-^f )r(i + 27) r(i + 7-^f )r(i-27) 
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Hence, as p — > 0, 



sin^(l + 27) \V^^^^[{P-R)p^-^i^ -{Q-S)p-^-^i^] ^ ^ ^ 



The physical solution ^ in Eg. ([58]) must belong to the domain of self- adjoint ness given by Vz{Hp). In fact behavior 
of the elements of the domain Vz{Hp) determines the boundary conditions for the system. If ^ ^^(^p), then as 
p — > the coefficients of r^"*^^/^^ and r~^~*^^/^^ in Eqs. ([5Q|) and (|58|) must match. Let us define 



m+3(P+ + i?+) = Xie'^^ 



and y m + ^(Q+ + 6'+) = X2e*^ 
Now comparing the terms of Eqs. (|5Q|) and (|58]) we obtain 

^2 \7p^^ xicos (01 



^^+m2y (5 + 5' X2COS ((/)2 ^ 2J 

Using the expressions of the constants P, Q, and S and the above equation we finally get 



d2 



- m 



(i + 7-+)r(i-7-+)r(i + 27) 



(59) 



= ^^^"^ + . (60) 

X2COS (02 + f ) 

Eq. (|6Q|) determines the spectrum in terms of the system parameters and the self-adjoint extension parameter z. Each 
choice of z corresponds to a different boundary condition described by the domain Vz{Hp) and leads to an inequivalent 
quantum theory. However the choice of z for a particular system is determined empirically as the theory cannot predict 
its value. For a special choice oi z = zi such that 02 + ^ = f , we have 

7- — = -P, P = l,2,3, (61) 

V 

This leads to the spectrum obtained in Eq. (|2T]) for < 7 < ^. For another special choice of z = ^2 such that 
^1 + f = f ' 

Ea 

-7 = -P, p = 1,2,3,.... (62) 

Though Eq.(|6Q|) cannot be solved analytically, from a typical plot of f{E) it can be obtained numerically. From Figure 
(3a) we can see when z changes from —0.8 to 4, the bound state energy changes from 0.915 to 0.92. Again when z 
changes from 4 to 0.1, the bound state energy changes from 0.92 to 0.942. Now calculating the contribution to LDOS 
coming from a single angular momentum channel j = f for these three different values of bound state energy, we 
observe in diagram 3(6) how the r dependence of LDOS varies with different values of z. In Figure (3c) we have shown 
how the bound state energy depends on the topology of the system for a particular angular momentum channel j = \ 
and three self-adjoint extension parameters z = —0.8,4 and 0.1. From the Figure (3c) we can see when n changes 
from 3 to 1, for the self-adjoint extension parameter 2: = 0.1 the bound state energy changes from 0.9935 to 0.9955. 
Calculating the contribution to LDOS coming from a single angular momentum channel j = \ for these two different 
values of bound state energy, we observe in diagram 3((i) how the topology of a system affects the r dependence of 
LDOS . 

In the scattering state sector where \E\ > \m\ and r] = —ik the general solution of Eqs. (fT4|) and (fT5]) that lead to 
physical scattering states are given by 



F{p) = Pi(fc)M (7 - 1 + 27, + Qi(fc)p-2w (_^ _ 1 _ 27, . 



(63) 
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FIG. 3: (a)Plot of f{E) is shown for system parameters i = f , n = 1, a = 1.48 and m = 1. The three horizontal hne 
corresponds to the three different values of the self adjoint extension parameter z = 4, 0.1, —0.8. (b) Dependence of LDOS in 
the bound state sector of the gapped graphene cone on the distance r from the external charge impurity is shown for three 
different values of bound state energy corresponding to three different values of self adjoint extension parameter. Here the 
contribution from angular momentum channel j = | is shown and the system parameters are n = 1, a = 1.48 and m = 1. 
We have assumed d=l. (c) Plot of f{E) is shown for two different values of n (1 and 3) with system parameters j = |, 
a — 0.29 and m — 1. The three horizontal line corresponds to the three different values of the self adjoint extension parameter 
z — 4, 0.1, —0.8. The solid lines correspond to n = 3 and the dotted lines correspond to n = 1. (d) Dependence of LDOS on 
the distance r from the external charge impurity is shown for two different values of bound state energy corresponding to two 
different values of n with the self-adjoint extension parameter z — 0.1. Here the contribution from angular momentum channel 
j — \ \s shown and the system parameters are a — 0.29 and m = 1. 



From Eq.([Tl]) we have 



G{p) 



lll^A(fc)Mfl + ,-^,l. 



)i(%-'^M 1-7 



Ea 



-,l-27,/0 



(64) 



Substituting these expressions of F and G in Eq (pTj) and Eq (p!2]) we get the upper and lower components of the 
wave function. Then from the asymptotic form of the wavefunction, identifying the incoming and outgoing waves, the 
scattering matrix and the phase shifts are obtained. Now, to find a relation between the constants Pi{k) and Qi{k), 
we consider the short distance limit of the wave function. The domain of self- adjoint ness of the Hamiltonian Hp is 
given by Vz{Hp) = V{Hp) {e*^^+ + e~*^^_}. In the limit r ^ 0, an element of the domain Vz{Hp) can be given 
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FIG. 4: (a)Phase shifts in the gapped graphene cone is shown for three different values of the self adjoint extension parameter 
z = 3, —0.5, —5 where the system parameters are n = 1, j = |,a = 1.48, and m = 1. (b) Scattering phase shifts are shown 
for different angles of the gapped graphene cone with the sae parameter z — —0.5 and system parameters j = ^, a = 0.29 and 
m—\. 



by 



5f = C(e'4«'+ + e-*t5f_), 



(65) 



where C is an arbitrary constant. 

After using the relation (|65]) and matching the coefficients of appropriate powers of r at both sides in (j65j), we get 
the following two conditions 



c 



sin7r(l + 27) 



m - 



d 



m 



(66) 



and 



(27?)- 



7 + 



»7 , 



Qi{k)yM + E 



-C 



(67) 



where z is the self-adjoint extension parameter. The equations Eq. (|66l) and Eq. (|67|) yield 



(^^ + ^+7-^)Pi(fc) 

{1^ + ^-7- ^)Qiik) 



m 



(Q+ + 5+) + , /m - ie-'^ (Q+ + 5+) 



-27 
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FIG. 5: (a) Dependence of LDOS on the distance r from the Coulomb impurity is shown for two different values of sae parameter 
z — —0.5,0.7 and a particular value of ^ = 4 and with j = |, n = 3, a — 0.29 and m — 1. (b) Effect of topology on r 
dependence of LDOS is shown for sae parameter z — —0.5, E — A and system parameters a — 0.29, n = 1,3 and m — 1 
considering contribution coming from the angular momentum channel j — \- (c)Energy dependence of LDOS is shown for two 
different values of sae parameter z = 5, — 0.8 at a distance r — 1 from the external Coulomb impurity. The system parameters 
used for the plot are a = 0.29 and m = 1 and contribution coming from the angular momentum channel j = | is considered, 
(d) Effect of topology on the energy dependence of LDOS is shown for sae parameter z = —0.8, angular momentum channel 
3 — \ and system parameters a — —0.29, n = 1, 3 and m — 1 



Xi cos( 



X2 COS(02 



-{2r^,f\2r^)-^\ (68) 



where we have defined y m + ^(^+ + R+) = Xi^ ^^^d ym-\- + S-^) = X2^^^'^ - Using the relations between 

the constants given before and (|68]), the scattering matrix can now be written as 



XI cos((/)i + f ) (-o^, \27/o^N-27 I+/2 r r(l+27) , n T{l-2-f) 
X2cos(02 + f)^^^l>' ^^^^ r(l+7-i^)^ ^ T{l-l-i^f 
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where 



iaE _ I iaE 

k £ — I ^ k (JVi) 

ima ' J'^ ima ' ^ / 

k ^ k 



The expression in (|69]) gives the S matrix for gapped graphene for the parameter range < 7 < ^ . For this range of 
7, the appropriate boundary conditions for which the Hamiltonian is self-adjoint and the corresponding time evolution 
is unitary, requires the introduction of an additional real self-adjoint extension parameter which labels the allowed 
boundary conditions. The phase shifts and the S matrix depend explicitly on z. For each value of z (mod 27r), we 
have an inequivalent set of the scattering data. Practically realizable value of z should be determined empirically as 
it cannot be determined analytically. 

In Figure 7(a) we have shown the energy dependence of scattering phase shifts for three different values of the 
self adjoint extension parameter z. We can see from the Figure that the region |£^| < \m\ is characterized by the 
sharp discontinuous oscillations which indicate the appearance of discrete bound states. In the other region where 
1^1 > |m|, scattering phase shifts behave qualitatively in the same manner for different values of z but they are clearly 
distinguishable from each other. In Figure 7{h) we have observed the effect of topology on the energy dependence of 
scattering phase shifts for a particular self adjoint extension parameter. It should be noted that the Figure shows the 
effect of topology considering only two values of n(n = 1,3) because the value of 7 remains in the region < 7 < | 
for only those two values. During the analysis we have always restricted the obtained results to the parameter range 
< 7 < ^ through the appropriate choice of system parameters. In the plots we have assumed that d = 1. 

In Figures 8(a) and 8(6) we have shown the dependence of LDOS on the distance r from the external Coulomb 
impurity placed at the apex of the gapped graphene cone where the energy is fixed at a value E = 4m. We have used 
the following expression given in equation Eq. ([7T]) for LDOS during the plotting. 

Mi5,r) = -A-^|*«(fc,r)|2. (71) 

3 

For the numerical calculation we have used the Equations (pTj) . (p!2]) . (|63|) . (|64|) and (|68|) . In Figures 8(c) and S{d) we 
have plotted the energy dependence of LDOS at a distance close to the charge impurity (r = 1). From these Figures 
we can observe that LDOS depend on the values of self adjoint extension parameter z and also on the topology of 
the system. Therefore measurement of LDOS using scanning tunneling microscopy can give us information about the 
self adjoint extension parameter and the topology of the system. 



IV. PROPERTIES OF GAPPED GRAPHENE CONE WITH SUPERCRITICAL COULOMB CHARGE 

In the supercritical region the radial part of the Dirac equation obeyed by the gapped graphene cone appears to 
be the same as Eq. (fT3|) . The only difference is that in this region 7 is always imaginary as the Coulomb potential 
strength a exceeds the value of u. We denote 7 = iA = \J o? — . Then from Eq. (fT3|) we have 

dF . ^ olE . ^ . ma.^ , 

p— + {t\ )F-{i^+ )G = 0, 72) 

dp T] 77 



and 



P— + - p + )G + {-ly + )F = 0. (73 

dp T] T] 



Substituting the expression of G from Eg. ([72]) in Eg. ([73]) we have 



cyE 

pF" + (1 + 2i\ - p)F' - {iX - —)F = 0. (74) 

V 

Solving the differential equations we can obtain the low energy eigenstates of the gapped graphene cone using two 
different boundary conditions. In the next section we shall consider a regularized Coulomb potential and obtain the 
quasibound state energy spectrum and then we shall observe how the nonzero mass and cutoff parameter affects the 
critical charge of the system. We shall repeat the same calculations with the zigzag edge boundary condition also. 
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A. Regularized Coulomb Potential 



In order to observe the supercritical effect of external Coulomb charge on the gapped graphene cone we shall first 
consider a regularized Coulomb potential because in that case we are allowed to extend the bound states until the 
negative continuum E = —m is reached ^581 [cij. The regularization of the Coulomb potential given by 



V{r) 



-a/r, r > a 
-a/a^ r < 



(75) 



where the Coulomb charge is placed at the apex of the gapped graphene cone, a is the minimum distance of the Dirac 
electron from the apex and it is of the order of the lattice parameter. The Dirac equation for gapped graphene cone 
is solved for these two different regions. 

Let us first consider the region r < a. In this region the Dirac equation is given by 



0. 



Eq. (|76p gives the following two coupled first order differential equations. 

}^{'\r) + - m + ^) P[^\r) = 



and 



E + m + ^j Q[^\r)=0. 



Substituting the expression of Qi\r) from Equation (|78|) in Equation (|77|) we have 



m > r 



= 0. 



A general solution of this Bessel equation is given by 



Using Eq. (|8Q|) we obtain 



(E+^-m) 
" y (E+^+m) 



(76) 



(77) 



(78) 



(79) 



(80) 



(81) 



Now we consider the region r > a. In this case the coupled 1st order differential equations and the second order 
differential equations obeyed by the Dirac fermions will be the same as Equations ([72]) . ([73]) and (f74|) . Then using 
the regularity condition at infinity we have [tE} 



and 



G 



= { 



F = U iX 



ma 



Ea 



1 + 2iX,p 



V V 



Ea 
V 



, 1 + 2zA,p 



Therefore the upper and lower components of the Dirac wave function will be given by 



(iA-i) 



Ea 



U (iX , 1 + 2iX, p + u]U (l^iX ,1 + 2iX, p 



Ea 



T] 



(82) 



(83) 



(84) 
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FIG. 6: (a) Bound state energy spectrum with regularized potential is shown. Here the blue line represents Arg[r(l + 2iX)\ 
and the dark red line represents RHS of Equation (|89|) ). (b) Dependence of |^(r)|^ on the distance r from the charge impurity 
placed at the apex of the gapped graphene cone is shown for a particular energy E — 0.96m obtained from the plot of the 
bound state energy spectrum. 



and 



U iX 



Ea , \ / ma \ / , Ea 

-,l + 2iA,p) - ( — -u]U(l^iX-—,l 



To determine the bound states we use the continuity condition of the wave function at r = a. It is 



-2i\p 
given by 



(85) 



Using this condition given in Eq. ()86p we have 

[/(iA- ^,l + 2iA,27?a) 



{^-^)u(l+i\-^A + 2i\, 2r,a) 



where 



When a 



f{E) 



l{m + E){E+^-m) Jl^+ll{^y{Ea + - m^a^) 
, (m - E){E + f + m) J^^_,^(^{Ea + a)^ - mW)' 

0, the continuity condition given in Eq. (|86]) will be satisfied when 

iX 



T{l-zX-^) 



u- ^{m-E)^iX 
V- ^{m-E) -iX 



iX 



r(i + 2iA) 
r(i - 2iX) 



OTf{E)=friE)=ATg[T{1^2iX)] 



(87) 



where p is a positive integer and 

Ea 



fr{E) = Arg 



F 1 + iA 



V 



■ Arg 



a 

V ym — E) ^ iX 

V 



Xln{2r]a) + Arg 



ly — a— — iX 



- pTT. 



Eq. (|89]) gives the bound state energy spectrum in presence of a regularized Coulomb potential for all possible 
situations. 
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Now we concentrate to the physically interesting cases where \E\ ^ m and the Coulomb potential is near to its 



critical value i.e a ^ ar 



for a particular n. Then from the Equation (|89|) we have 



In 



^[J\u-^\{^)-J\u+i\{^)] 



■ ^{1 



Here 



r{x) 



and 77 



— m^. For ^ << 1 we have up to the terms of order 



En 



2 En 



2a 



V-il 



(^) 



(1-0 

2iy 
(1 - 



T 



1+1 



(91) 



Thus we can see the effect of the nonzero mass on the bound state energy spectrum. 

The mass affects the critical charge of the system. For a regularized Coulomb potential the bound states can dive 
into negative energies. Here the critical charge refers to that value of Coulomb potential for which E = —m. Then 
for the region near critical potential we have 



2v\og^ [2mvCa] 



where 



C = exp 



(92) 



(93) 



From Equation (|92]) we can see when ma 0,ac ^ ^ which agrees with the result obtained for massless case [76|. 
The dependence of critical charge on the nonzero mass and cutoff parameter are shown in Fig.® for different opening 
angles of the gapped graphene cone. From the Fig.© it is clear that the topology affects the critical charge of the 
system and the nature of their dependence on the product of mass and cutoff parameter of the system remains almost 
same. 



B. Zigzag edge boundary condition 

To find out the energies of the stationary states formed in the supercritical region now we use the zigzag edge 
boundary condition ^"^(a) = 0, where a is a distance from the apex, of the order of the lattice scale in graphene. In 
order to proceed we shall first solve EQ.([7i|) and obtain two linearly independent solutions Fi{p) and F2{p) which are 
regular at p = 0. They are given by [75|] 

Ea 

Flip) =AiM(iX ,l + 2iA,/)) (94) 



and 



Then from Eg. ([72]) we have 



and 



F2{p) = A2p-^'^M{-iX - — , 1 - 2iA, p). (95) 

V 



(iX-—) Ea 
Giip) = \ Z^ ArM{l + , 1 + 2iA, p) (96) 



G2{p) = -^^^^^Mp-^'^M{1 -iX-^,1- 2i\, p). (97) 
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FIG. 7: (a)Bound state energy spectrum with zigzag edge boundary condition is shown. Here the blue hne represents Arg[r(l + 
2iX)] and the red hne represents RHS of Equation ()100p ). (b)Dependence of |^(r)|^ on the distance r from the charge impurity 
placed at the apex of the gapped graphene cone is shown for a particular energy E — 0.94m obtained from the plot of the 
bound state energy spectrum. 



The solution satisfying the zigzag edge boundary condition can be given by 



^^^(r) = vW^e-«p^-2 [{Fi(a) - Gi(a)}{F2(r) - G2(r)} - {Fi(r) - Gi(r)}{F2(a) - G2(a)}] 



(98) 



The square integrability condition of the wave function indicates that as p — > cx3 the diverging part of the wave 
function must vanish. Therefore we have 



/ P { —M^ — . 



r(l - 2iA) 



(iA-f^)' 



r(i + 2iA) 
r(iA-f^) 



HA-f^) 
- (p-i- 



-2iAln(2i7a) 



2Ea \ 



This gives the condition 

r(l - 2iA) 



THA-f^) \ 

r(iA 



(iA-f^)' 



or j{E) 



Eg 



Eg -. 



„2iAln(277a) 



r(-iA 

or = = Arg[r(l + 2iA)] 



r(i + 2iA) ( 

'v+^{m + E)-i\' 



HA-f^) 



-2iAln(2r7a) 



0. 



(99) 



r(l + 2zA) 

r(i - 2iX) 



(100) 



where p is a positive integer and 

Arg 





+ Arg 















a 
-( 



iX 



Xln{2r]a) -\-p7r. 



This is the condition directly obtained from our analysis and obeyed by the bound state energy spectrum for all 
possible situations. 

We can compare this bound state energy spectrum for zigzag edge boundary with the spectrum obtained from the 
regularized potential case and observe how the two different boundary conditions affect the spectrum. From Fig.® 
we can observe that the bound state energy spectra and the probability amplitude of the wavefunction gets affected 
by the boundary conditions though their nature remains same. 

Now we again concentrate to some physically interesting cases with some approximations. First we consider the 



case where \E\ m and the Coulomb potential is near to its critical value i.e a 



an 



^mm 



for a particular n. 
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FIG. 8: (a) Bound state energy spectrum with zigzag edge boundary condition and regularized Coulomb potential. Here the 
blue line represents Arg[r(l + 2iX)] and the dashed and the solid line represents RHS of Equations (|1QQ|) and (|89l) respectively, 
(b) Dependence of |^(r)|^ on the distance r from the charge impurity placed at the apex of the gapped graphene cone is shown 
for both the zigzag edge boundary condition and regularized Coulomb potential. The values of energy are obtained from the 
bound state energy spectrum. From Fig. (a) we can see that for zigzag edge boundary condition a possible bound state energy 
is ^ = 0.9895m and for regularized potential a possible energy is ^ = 0.995m. Here we have given the probability amplitude 
plots for these two particular energies. 



Then from the Equation (jlOOp we have 



In 



-2iEr, 



m 



-.a] =2^(1)-^ 



^2 



T' 



(101) 



Like the previous section here also il){x) 



mi 

' 



and T] = —i^l — m? . For ^ << 1 we have up to the terms of order 



En 



or Er, 



2 En 



2 En 



2a 



2ip{l) -ip{-iu) 



(1 - i) _ pn 
2v A 

(1 — i) pn 



^ m 






}-e; 







4u 



2y 



A 



1 



m f 1 \ m? \ iv 



(102) 



Thus we can see the effect of the nonzero mass on the bound state energy spectrum. 

Like the critical charge obtained for regularized Coulomb potential, here also with zigzag edge boundary condition 
we can see that the mass affects the critical charge of the system. Proceeding as before in this case we have 



TT 



2v\og^ [2mz/aexp ( - 2^/^ ( 1 ) ) ] 



(103) 



Here also from Equation (jlQ3p we can see when ma ^ 0, ac ~ which agrees with the result obtained for massless 
case |76]. The dependence of critical charge on the nonzero mass and cutoff parameter has been shown for different 
opening angles of the gapped graphene cone. Also we have compared the dependence for two different boundary 
conditions. 

From Fig.(|9]) we can see that for a zigzag edge boundary condition the critical charge vary with ma more rapidly 
than it varies for a regularized Coulomb potential. Thus the two different boundary conditions affect the dependence 
of critical charge of a gapped graphene cone on ma. 
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FIG. 9: Dependence of critical charge on the nonzero mass and cutoff parameter are shown for both zigzag edge boundary 
condition and regularized Coulomb potential for different opening angles of the gapped graphene cone. The dotted lines show 
the dependence for zigzag edge boundary condition and the solid lines show the dependence for regularized Coulomb potential 

V. CONCLUSION 

In this paper we have described the low energy dynamics of massive Dirac fermions in a gapped graphene cone in 
the presence of an external Coulomb charge impurity. The graphene cone can be equivalently described as a graphene 
plane together with a flux tube whose gauge potential is chosen to produce the required holonomies. The given system 
thus consists of a gapped graphene plane together with a combination of a Coulomb charge impurity and a flux tube 
passing through it. The strength of the Coulomb charge can be sub or supercritical. 

The combination of this topological defect as well as the charge impurity results in short distance interactions, the 
effect of which cannot be incorporated as dynamical terms in the Dirac equation, valid in the low energy limit. For a 
sub critical charge impurity, we show that the effect of these interactions can be modelled through appropriate choice 
of boundary conditions, which are determined by imposing the requirement of a unitary time evolution. While there 
is a very large class of allowed boundary conditions, it turns out that they can be labelled by a single real parameter. 
It is this parameter through which the effect of the short range interactions enter in the analysis presented here. This 
parameter cannot be determined from theory alone. However, we have shown that the observables such as LDOS, 
scattering phase shifts and bound state energies depend explicitly on this parameter. As mentioned before, a similar 
situation arose in the context of Dirac fermions in a plane in the presence of a cosmic string. These two situations are 
not identical, but similar quantum subtleties arise in both contexts and the case of the gapped graphene considered 
here is more amenable to empirical analysis. 

The supercritical regime of the external Coulomb charge is characterized by quantum instabilities. Here we have 
analyzed the effect of topological defects in the supercritical dynamics of gapped graphene. This problem has been 
analyzed with a regularized Coulomb potential as well as with a zigzag boundary condition. We have shown that the 
quasibound state spectra and the probability amplitude depend exphcitly on the number of sectors removed from a 
planar graphene to form the cone. In addition, the running of the critical charge as a function of the product of the 
Dirac mass m and the cutoff parameter a has been obtained. Though the nature of the dependence is similar for 
both the regularized Coulomb potential and the zigzag boundary condition but in the latter case the critical charge 
increases more rapidly with ma than the former case. 

It would be interesting to study the effect of sample topology on the electron-electron interactions in graphene and 
the associated gap equation, which is currently under investigation. 
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